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Abstract 

Cosmological gamma ray bursts are very likely powerful sources of high 
energy neutrinos and gravitational waves. The aim of this paper is to re- 
view and update the current predictions about the intensity of emission in 
this two forms to be expected from GRB's. In particular a revised calcula- 
tion of the neutrino emission by photohadronic interaction at the internal 
shock is obtained by numerical integration, including both the resonant 
and the hadronization channels. The detectability of gravitational waves 
from individual bursts could be difficult for presently planned detectors if 
the GRB's are beamed, but it is possible, as we have proposed in a paper 
two years ago, that the incoherent superimposition of small amplitude 
pulse trains of GW's impinging on the detector, could be detected as an 
excess of noise in the full VIRGO detector, integrating over a time of the 
order of one year. 
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1 Introduction 



The Gamma Ray Bursts (GRB) are intense flashes of soft 7-rays, serendipitously 
discovered by the Vela test ban treaty (Klebesadel et al. 1973) and immediately 
confirmed by the Soviet Konus satellites (Mazets et al. 1974). The phenomenol- 
ogy of GRB's has been discussed in several reviews (see e.g. Djorgovski et al. 

(2001) and references therein). The most important breakthrough in the GRB's 
observations has been the discovery in 1997 of a long lasting X-ray "afterglow" 
of the burst GRB970228 by the Beppo-Sax Italian-Dutch satellite (Costa et 
al. 1997), which has allowed the identification of its host galaxy (van Paradijs 
et al. 1997, Sahu et al. 1997), making the study of GRB's a "multi- wavelengths 
adventure" (Klose 2000). 

According to a widely accepted (and acceptable) interpretation, the prompt 
7-rays burst and the delayed afterglow, are both powered by the dissipation 
of the kinetic energy of a relativistically expanding "fireball" (or jet), whose 
primal cause is not known. This mechanism, qualitatively suggested by Cavallo 
& Rees (1978), has been refined by Paczynski (1986) and put in its final form of a 
standard model by Rees & Meszaros (1994) (for a recent review see e.g. Meszaros 

(2002) and references therein) . The obvious expectation is that such a powerful 
highly relativistic explosion should be accompanied by emission of neutrinos 
and gravitational waves, thus becoming the herald of the new astronomy. 

The aim of this paper is to review and update the current predictions about 
the intensity of emission in this two forms to be expected from GRB's. In the 
following I will first review in ^ and 13 the essential quantitative information 
on the energetics and the large scale distribution of the sources, which has been 
obtained from prompt and afterglows observations. In 21l present a new calcu- 
lation of the photohadronic neutrino emissivity in the fireball/jet system from 
which predictions about the observability of the diffuse neutrino background 
can be derived. In f|Sll reexamine the possibility that a similar stochastic back- 
ground could be observed as an additional noise component in the long term 
integration of the signal from gravitational wave detectors. 

2 Energetics 

The total energetics involved in the GRB's is the most important parameter for 
my subsequent emissivity estimate. A direct estimate of the rest frame energy 
output can be obtained only for bursts with known redshift. According to a com- 
pilation by J. Greiner, available at http : / / www . mpe . mpg . de/~ j eg/ grbrsh . html, 
the number of GRB with reliable redshift measurement of the host galaxy red- 
shift, is about 20. 

The intrinsic isotropic bolometric (from 0.1 keV to 10 MeV) energy release, 
estimated from the fiuence of 17 burst in the 30 keV -2 MeV range, spans from 
5.6 X 10^1 erg (GRB990712) to 2 x 10^^^ erg (GRB990123), wfth a median of 
= 1.8 X 10^^ erg and a dispersion of 0.8 dex(r.m.s.) (Bloom et al. 2001). 
However, it is to be emphasized that this value applies only to GRBs with 
measured redshift. Therefore observational biases in burst detection and redshift 
determination, might obscure the true underlying average energy output of the 
typical GRB. However, keeping in mind this caveat, one can assume that the 
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kinetic energy input to the expanding relativistic fireball is of the order of 



where r}^ is the fraction of the fireball kinetic energy that is converted into 7rays 
(Guetta et al. 2001). 

The afterglow light curve of many bursts shows an achromatic break, cur- 
rently interpreted as the evidence for a beamed emission (Sari et al. 1999, Rhoads 
1999). In this case the emitted 7-ray energy S-!^^^ is obviously reduced by the 
factor (Frail et al. 2001) Si^*^* = fjet^^"^" where fjet = 1 — cos jet- In a recent 
paper (Bloom et al. 2003) present a complete sample of 29 gamma-ray bursts 
(GRBs) for which it has been possible to determine temporal breaks (or lim- 
its) from their afterglow light curves. Incorporating realistic estimates of the 
ambient density and propagating error estimates on the measured quantities, 
in agreement with the previous analysis of a smaller sample, the derived jet 
opening angles of those 29 bursts result in a narrow clustering of geometrically 
corrected gamma-ray energies about (S^^*) = 1.33 x 10^^ ergs, with a burst-to- 
burst variance of 0.34 dcx. 

Three consequences can be drawn, if this is true: 

1. The central engines of GRBs release energies of the order of 

£0 « few X 10-2 (^^^ Mqc^ 

This value is indeed comparable to the kinetic energy output of ordinary 
SNe, suggesting a connection between the two phenomena, also supported 
by the observation of bumps in the afterglows light curves (Bloom 2003); 

2. The large spread in observed fluence and peak luminosity of GRBs is 
mostly due to a spread in the opening angle of the beam; 

3. Only a small fraction of GRB's are visible to a given observer and the true 
GRB's rate is 

^GRB — ^GRB ^ {fjet ) 

In practice from Table 2 of the paper by Bloom et al. (2003) one can 
guess that the true GRB's rate should be about 300 times larger than the 
observed one. 

It is clear that such a radical descoping of the energetics associated with 
the GRB's explosion, has a fundamental impact on the physical interpretation, 
but has little or no impact, at the first order, in the predictions of their high 
energy neutrino emission. In fact we expect that the aperture of the beam of 
neutrinos emitted by photohadronic processes will be the same of the beam of 
7-rays emitted by electrons. 

On the contrary the gravitational wave emission by GRB, if present, will be 
in any case isotropic. We observe that the beaming could reduce the probability 
of detecting individual events, because the amount of energy to be radiated is 
strongly reduced. But this should have no impact on the prediction of a possible 
cosmological background originated by the GRB, because in case of beaming we 
should assume that the rate of events is increased by the same factor. Therefore 
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we have that the total energy released is in the case of beaming will be in the 
average 

^GW ^ ^GRB — ^GW ^ ^GRB 

As we said already the above given estimate can be applied, only to the 
bursts with observed afterglow. Since afterglows have been observed only for 
long lasting bursts with duration 190% > 2 s, it could be argued that the short 
bursts could be powered by a totally different central engine, whose energetic 
could be on a lower scale. Some indications in this direction is given by the fact 
that the observed fluence of the short bursts is significantly lower then the one of 
long bursts (Tavani 1998). However one could argue back that the difference in 
the average fluence could be due to a markedly different large scale distribution 
of the two classes of bursts. In the absence of a any direct distance estimator, 
only statistical indicators of the spatial distribution of GRB can be used to test 
this hypothesis. A recent analysis of the 4B burst catalog (Schmidt 2001) shows 
that the short bursts have essentially the same characteristic peak luminosity 
of the long ones, but their local space density is around three times lower. This 
applies to bursts shorter then 0.25 s, therefore the total energy output of the 
short bursts should be 1-2 order of magnitude lower then the one of the long 
burst. 



3 Large scale distribution 

The common wisdom about the large scale distribution of the bursts is that, 
being a phenomenon originated by the collapse of short lived massive stars, 
they should track closely the Star Formation Rate (SFR) cosmological distri- 
bution. In a recent paper Porciani & Madau (2000) have used three different 
SFR derived from other astronomical data to fit the observed log N-Log P 
distribution for the bursts. Their results imply about 1-2 GRB's per million 
SN type II and a characteristic isotropic luminosity of the bursts in the range 
3 — 20 X 10^^ ergss~l. The best fitted models of the rate per unit comoving 
volume are shown in Fig. ^ 

The luminosity-variability relation(Fenimore & Ramirez-Ruiz 2000) is a new 
entirely empirical connection, which has been proposed in analogy to the well 
understood luminosity-period relationship of the cepheids. Having calibrated 
over a sample of 8 GRB's with known redshift, the authors have applied this 
fit to derive the redshifts for 220 bursts. Applying non-parametric statistical 
techniques to the 220 Gamma-Ray Burst (GRB) redshifts and luminosities de- 
rived in this way, it has been proposed (Lloyd-Ronning et al. 2001) that there 
exists a significant correlation between the GRB intrinsic luminosity and red- 
shift, which can be parameterized as (17-y/47r) L-y cx (1 -f z)^-^^°-^, where z is the 
burst redshift. In addition, this analysis supports a co-moving rate density of 
GRBs that continues to increase for [1 + z) > 10. We have plotted in the Fig.Q] 
the different possible large scale distribution of the GRB's, which is compatible 
with the observed log N-log P distribution. 
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Figure 1: 



4 High energy neutrinos 

In the general framework of the relativistic firebaU model of GRB's, the prompt 
7-ray emission is originated by synchrotron emission of relativistic electrons 
accelerated at the internal shocks, while the afterglow is due to the interaction 
of the external shock with the ambient material (Piran 2000). The radius of 
the inner shock (or shocks in case of multi-peaked emission) should be of the 
order of i?s ~ 2 c At, where At is the smallest timescale observed in the burst 
emission, usually in the range 1-10 ms. 

A lower limit to the bulk Lorentz factor of the expanding fireball can be de- 
rived from the 77 — > e"*" e~ self absorption in the burst at the highest observed 



7-ray energy 
2 rrip c 



100 MeV. In fact given the threshold of the reaction 

100. A more detailed 
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- 1 MeV we must have F^ > e™"^/V27H7? 
calculation (Lithwick & Sari 2001) shows that the optical depth 
where a is the low energy spectral index of the prompt 7-ray emissioin. There- 
fore a bulk Lorentz factor F^ > 300 is required for the emission of ~ 100 MeV 
7-rays. 

It has been shown (Vietri 1995, Waxman 1995) that the high magnetic fields 
and relativistic velocities of the internal shocks, can accelerate protons up to 
10^^ — 10^° eV, corresponding to the range of the UHE cosmic rays. From 
first principles one can derive the rate of acceleration from the relativistic shock 
£p/ep ~ ;^_BF^, thus even in a very short time, smaller than the coherence 
time of the 7-ray emission, a huge maximum energy 
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can be reached. 

However the shock acceleration mechanism is efficient only when the radius 
of the shock is larger then the Larmor radius ~ cep/eB of the accelerated 
particle. Therefore the more stringent limit will be ^ Rs, that assuming 
Rs ~ 2cFjAi gives e™"^ < eSF^At. But even in this case one can get 
gmax ^ ^q19 gy Tff[^]^ Held of the order 1 T. The typical fields required in the 
burst for electron synchrotron emission peaked in ^ 100 keV range are S> 1 T 
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(Lloyd & Petrosian 2000). Therefore the maximum energy imparted to the pro- 
ton wiU be limited, in practice, by proton synchrotron radiation and inelastic 
scattering over the radiation field. The latter interaction will copiously produce 
pions, that subsequently decay into neutrinos and TeV 7-rays. The pion yield 
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Figure 2: 

for the reaction pj ^ X at the internal shock has been calculated, under 
simplifying assumptions, by several authors (see e.g. Guetta et al. (2003) and 
references therein). All this previous calculations have included only the photo- 
hadronic production near threshold, where the photon scatters quasi-elastically 
with the nucleon that is pushed into an excited state corresponding to iV^ and 
A"*" baryonic resonances: 

^^^{ A+ }^{ 

These resonances are produced with high cross sections (~ 100 — 500 /xbarn) 
(Hagiwara et al. 2002). In addition to these resonant channels, pions will be also 
produced at intermediate energies by the hadronization of the photon 7 qq. 
The reaction proceed in practice as 

{qq)p^M,TT^,TT"r-- (2) 

The importance of considering both processes, and not only the resonant chan- 
nels, in photohadronic pion production calculations has been already stressed 
in the literature (Mucke et al. 1999). In addition to the arguments given in that 
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paper, one should consider that the photons which interact via the resonant 
channel should have energies in the range 



P ^ ^ max (ry\ 

In Fig. 121 we have plotted the cross section versus the energies, in the observer's 
frame, of the photon and the 7-ray. If the low-energy spectral index is a < 0, 
protons with > 10^'^ eV will produce baryonic resonances scattering on pho- 
tons with < 1 keV. Therefore we expect that the resonant photoproduction 
at high energies will depend strongly on the shape of the 7-ray spectrum, in a 
range well below the observed one. 

The theoretical expectation from the electron synchrotron emission (Sari & 
Piran 1999) is that the spectrum will level off due to electron cooling at energies 
£7 ^ 100 eV and be strongly suppressed by self-absorption for < 10"'^ eV. 
Following this indication we have assumed a photon energy distribution given 
by the form empirically proposed by Band et al. (1993), modified to take into 
account the cooling: 

£;a g(/3-a)SJ£;, < 

(e^) = A <( e" ^{li-o^y.lE, if < 

in the range 10""^ eV < < 100 MeV and zero outside this range. The constant 
A has been chosen in order to normalize to 1 the function ipi^-y) in the above 
given range. Under these assumptions the photon density at the source will be 



At 



(6,) (5) 



where Rs ~ 2c Fg At. 

The full expression for the m.f.p. of a proton with energy ep in a radiation 
field with photon energies e™*" < £7 < e™"^ in observer's frame is given by 



where, cjp-y is given by the sum of Breit-Wigner resonances and hadronization 
cross sections 

k=N+.A+ Ws-Mk) + 

, we can calculate the average number of inelastic interactions per proton 

I , AR _ L^\-^ 
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The probability of having one or more proton-photon inelastic scattering is 
= 1 — exp(77^) which is Pn- « t],^ only when 77^ <C 1. This is not the case, 
specially in the hadronization channel, because the inelasticity of the interaction 
is of the order of kinei — 40%. Therefore one has to take into account that the 
proton actually starts an small hadronic cascade in the considered layer. We 
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Figure 3: Fraction of energy lost by protons to pions, assuming a photon energy 
distribution given by Eq. |0J,with low-energy spectral index a = —1, high- 
energy index (3 — —2.25, a break energy Ei, = 511 keV, and a cooling energy 
Ec = 100 eV. For the fireball's parameters it is assumed = 10^^ erg s~^, Fb — 
300 and At = 10 ms. 

have integrated numerically Eq. © including all the resonant channels and the 
continuum hadronization cross section, approximated by a polynomial logarith- 
mic fit (Cudell et al. 2000). We have also approximately taken into account 
the possibility of multiple interactions, simply summing over the histories. In 
practice we have considered that the pions of a given energy e^r could have been 
produced either by a proton of energy Sp or by a proton which suffered n pre- 
vious scattering, with n < log(ep/e7r)/ log(l — kind)- Thus we have in practice: 

"ma x 

= E r^""-(re.) {l-exp[,7^(re.)]} (9) 

n=l 

where is ^ = 1/(1 — fci„ei), 7 the spectral index of the proton spectrum, and n^r is 
either the average multiplicity of the multiple pion production for the hadroniza- 
tion channel or the branching ratios in the case of the resonant channel. In Fig. 
|Slwe report the fraction of energy lost by protons to pions = 77^ x {tTi/tp) for 
comparison with other calculations. As can be seen the result is in reasonable 
agreement with the first estimate given by Waxman & Bahcall (1997) which is 
shown as an heavy dot-dashed line in the same figure. In Fig. 01 we report the 
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Figure 4: Typical neutrino spectrum for a burst with Ep = 2.2 x 10^^ erg with 
spectrum cx for 10^° eV < ep < 10'^" eV, exploding at z = 1 in a Universe 
with Hq = 71 kms^^Mpc^^, JIm — 0.3 and JIa — 0.7 and w — —1 (see 
text). The fireball parameters are the same of Fig. 13 This hypothetical burst 
(assuming = £p and — 0.04) would have a fluence at Earth in the 20-2000 
keV band of 3 x 10^^ ergcni"-^. 



energy distribution for the neutrino produced in the decay chain 

TT ^ and /i — > e j/g (10) 

In both cases it is important to take into account the pion (and muon) en- 
ergy loss, before decaying. This has been done solving analytically the coupled 
equations for decay and energy loss, which gives as a result when the dominant 
energy loss is the synchrotron emission: 

where e*^ = (mj — rn^)/2m7r is the energy of the neutrino in the muon rest 
frame, and e^rc is the energy of the pion for which the rate of energy loss is 
equal to the decay time. For e^r ^ Cttc Eq. ((TT|l gives the expected behavior 
(X (eTr/circ)"^, while for « e^c a pile up peak is generated. We can derive 
the spectrum of the muons from the same formula with the only difference 
that the integration limits arc in this case 2 6™!+^*) — — 2e™(i-'/3*) '^here 

e* = \J''^fi + ^Jf, ^ is the energy of the muon in the pion rest frame and /3* its 
velocity. In the muon decay the energy loss before decay will be more severe 
because the decay constant of the muon is longer (implying e^c << Ettc), but 
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Figure 5: Neutrino spectrum for different bulk Lorentz factor keeping all the 
other parameters constant. 

we can use a formula equivalent to Eq. For the sake of simplicity we 

have calculated the combined neutrinos and anti-neutrinos flux, because in the 
multiple pion production it can be assumed that positive and negative pions 
are equally produced. The marked difference between the spectrum of muon 
neutrinos and the one of the electronic neutrinos is due to the fact that the first 
type of neutrinos are directly produced in the pion's decay, while the electronic 
ones are produced only in the muon's decay, which have suffered a larger energy 
loss. In Fig. 0] we show for comparison the neutrino flux calculated by Guetta 
et al. (2003), which is also in good agreement with this result. In Fig. we 
report our prediction for the expected neutrino flux from the photohadronic 
channel, for different values of the bulk Lorentz factor. It is to be remarked 
that I do not find in the final flux estimate the strong dependance cx F^* 
which is argued by many authors (see e.g. Guetta et al. (2003)). The reason 
is that even if the average number of inelastic collision Ty^r calculated by Eq. 

is cx F^"*, the pion yield is 77^ cx Ty^r only when ?77r <C 1. As we said above 
when it becomes ?77r ^ 1 one should take into account the multiple regenerative 
interactions, which take place at decreasing energies, as has been taken into 
account using Eq. l(TT|) . In Fig. we report the diffuse cosmological background 
flux of neutrinos obtained from our emissivity estimate and the different large 
scale distribution models discussed in fj^l The neutrino fluence at Earth of a 
burst with redshift z has been estimated as 



where dh{z) is the cosmological luminosity distance, £p the total accelerated 
proton energy with spectral index —2, and Vj^ is given by Eq. Hll|l . The diffuse 
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Figure 6: Diffuse neutrino flux for the different models of GRB's large scale 
distributions shown in Fig.^ for a cosniological model with the same parameters 
of Fig. 0] The normalization of the curves has been obtained assuming an 
average total energy for accelerated protons in the burst rest frame of £p = 
10^^ erg. The band shown in this figure corresponds to ±0.8 dex variation of 
this average. 



flux of neutrinos has been calculated integrating 

oo 

^, X [ Rgrb[z) ^dV 

= J -jY^—yT,y[z,ey)—dz (13) 







where 



(1 + ^) 
dV _ cdl{z) 



dz {l + z)^H{z) 
being for Friedmann-Robcrtson- Walker metric (Hagiwara et al. 2002) 



(14) 



H 



(z) = Ho ^JnM{l + z)-^ + nK{l + zf + fiy(l + z)-3(l+-) (15) 



with fix = 1 ^ ^^_A/ — ■ 

Finally in Fig. we report our prediction for the fluence of the 7-rays with 
> 100 GeV produced by neutral pion decays. This mechanism has been 
discussed in the literature (for a recent review see Fragile et al. (2002) and 
references therein) as a possible source for VHE 7-rays potentially detectable by 
ground based experiments. It is well known that only one statistically significant 
detection has been reported until now by Atkins et al. (2003) using the Milagrito 
detector. This type of signal is expected to be detectable only for very energetic 
and nearby GRB's for two reason: 1) the strong self absorption due to 77 
e^e~ interactions with the soft 7-ray radiation field in the source, and 2) the 
absorption of the VHE photons by the cosmic microwave background photon 
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Figure 7: The predicted VHE 7-ray fluence of a very nearby bursts {z = 0.01). 

due to the same process. The first process is responsible for the suppression of 
the emission below « 1 — 10 TeV, while the second originates the cut-off at 1 



5 Gravitational waves 

Prompt gravitational waves emission has been predicted from GRB since 1993 
by Kochanek & Piran (1993) at a level to be detected by upcoming gravitational 
wave-experiments. As we have discussed in |2 the estimate of the kinetic energy 
required to power a beamed 7-ray burst is of the order of £"0 > 0.01 Mq c^, much 
larger then the total nuclear binding energy of few-Mo's stars. Therefore in 
any case the energetic of the bursts is compatible only with the collapse of a 
several solar masses object into a black hole (Woosley 1993). Any reasonable 
assumption regarding angular momentum then leads to a massive accretion disk 
that flows, on a viscous time scale, into a black hole that very rapidly approaches 
the Kerr limit. The dissipation the rotational and gravitational energy within 
the geometry of a thick torus is very likely to lead to jets, possibly by the 
purely electromagnetic mechanism proposed by Blandford & Znajek (1977). 
The energy irradiated in this case as gravitational waves (Stark & Piran 1985) 



where m is the mass of the torus, Mbh is the mass of the black hole, and e is 
the ellipticity of the torus. It is worth noticing that if the jets are accelerated 
by Poynting flux (Drenkhahn 2001) the amount of kinetic energy going into 
the burst £0 and the energy radiated in gravitational waves are strictly related. 
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because the GW luminosity will be Lqw oc {Qij Q*-'), being Q the quadrupole 
moment of the mass distribution of the torus, while the Poynting flux will be 

Le.m. oc {Mij M*-' ) being in this case M the quadrupole moment of the magnetic 
field embedded in the torus. Therefore we assume that the energy radiated 
as gravitational radiation is Saw = Vgw £o where rjow 1 is an unknown 
parameter. 

The energy flux (viz. energy per unit surface and unit time) of GW produced 
by a burst of intrinsic luminosity (in gravitational waves) Lqw at a red shift z 
is by definition ^gw = Lqw I '^t^ d\{z) if is the luminosity distance. In order 
to convert this flux into an adimensional amplitude we can use the classical 
formula (Shapiro & Teukolsky 1983): 

^ow^^{h\ + j^.) (17) 

where the average is taken over several wavelengths. The amplitude of the signal 
produced depends from the direction and the beam pattern of the detector. In 
the best case we have integrating over time and applying the Parseval' s theorem 

4-00 

2?;2/ N , IQttG {\ + z) Egw 
" ^H^"= ,3 4.rdi(z) 

— oo 

where h'^{uj) will be the Rayleigh power of the signal as a function of the fre- 
quency Lu. In order to estimate the order of magnitude of the amplitude of the 
GW signal we do not need a detailed shape of the spectral power density of the 
signal, but only the knowledge of the firsts two moments of the distribution, 
viz. uj and Au. In fact we can recast Eq. H18(l in the form 

+00 

[lu + Auj ) f h [uj)duj ^ — ^ — T^-^- — 

^ 'J c3 47rd2(z) 

— oo 

It is remarkable that rather natural physical assumption on the first and second 
moment of the unknown Rayleigh power distribution h?{uj) can be made. In 
fact we can assume that the first moment will be twice the keplerian angular 
velocity of the marginally stable orbit around the BH, namely 



where rs will be the Schwarzschild radius of the collapsed object, which is 
function only of the black hole mass, and ^ < /3 < 3 depends on the angular 
momentum of the BH, being minimal for maximally rotating BH, and maximal 
for non-rotating ones. 

The second moment will be the r.m.s. bandwidth of a wave packet that 
can be estimated Alo k, 2tt/ At where At is the timescale of the emission in the 
comoving frame. As we have seen in ^the 7-ray emission has a variability over a 
timescale At ~ 10 ms, therefore it is rather natural to assume that the timescale 
of the GW emission by the torus will be of the same order of magnitude. In 
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this case we have, if c Ai 3> rs, f^or the peak amphtude 



dL (z) hP'^-'^ « J Mbh VEawAt (20) 

The expected amphtude for typical values of beamed bursts is 

compatible with our previous estimate (Auriemma 2001). The predicted signal 
by Coward et al. (2002) for a "gentle chirp" spectral model, is well in agreement 
with the above value, when scaled by a factor 10^ to take into account the 
different assumption on the distance and the total energy, which are in that 
paper respectively cIl = 100 Mpc and Sqw — 3.75 x 10^'^ erg. The current 
upper limit for the search for emission of GW in correlations with GRB with 
bar detectors (Tricarico et al. 2001, Astone et al. 2002) has given the upper 
limit hjfujs < 1-5 x 10^^*. It is clear from this estimate that the probability 
of detecting a single burst, even under optimal is very low, unless rjcw 2> 
1. However even if the individual burst could not be detected, it is to be 
remarked that in case of beaming the rate of explosion could be very large 
(« 1500 per day), therefore the stochastic accumulation of signal integrated 
over a long time could emerge from the noise, as we have shown in our previous 
paper (Auriemma 2001). The same type of calculation, but with more optimistic 
assumptions on the energetic of the source, as we said before, has been repeated 
by Coward et al. (2002). The energy flux of GW produced at Earth from 
a cosmological distribution of sources can be calculated integrating over the 
cosmological distribution, as we have done for the neutrino in the Eq. (|12|l . 
namely 



(1 + z) A'Kd\{z) dz 



-^^ '^LZ^T ^-^dz (21) 



where dV/dz is given by Eq. ifTljl . Applying again Eq. ((T7|l we have: 

hl) = ^_T^^'4^dz (22) 



16nG \ + At: J {l + zfH{z) 



On the L.H.S. of this equation we have the amplitude of the wave that invests 
at a certain instant the detector. We have seen in the previous section that 
each of this burst will not have a detectable intensity, but if we average over an 
observation time T long compared to the GW burst duration but short compared 
to Hubble time scale (typically one year) we have a signal 

-oo 

that will be detectable if the power spectral density is greater then the power 
spectral of the noise, averaged over the same observation time. The uncorrelated 
superimposition of bursts of gravitational waves will be well approximated, for 
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Figure 8: Stochastic signal integrated over one year from cosmological GRB 
{^GW — 0.01) with large scale distributions of Fig.^ multiphed by an average 
beaming factor (/~j) = 300. The dotted line is the noise estimated in the 
VIRGO experiment (see text). 



the central limit theorem, by the superimposition of redshifted gaussian distri- 
butions. Therefore the power spectral density of the signal can be estimated by 
the integral 

C2w2 7 72^(7^ (1 + Z)2 77(z) 



where ujq will the first momentum of the distribution of the frequencies uj(Mbh) 

given by Eq. ((11^ . namely luq — Jj^.j^ u>{Mbh) i^{MBH) dMsH where ^{Mbh) 
is the normalized mass distribution function for the black holes, and the 
second momentum. In the lack of any better guess we will assume that ^{Mbh) 
is flat for 4M0 < Mbh < 14 M©, thus having luq = 27r 10^ rads^^ and = 
27r200 rads-^ 

Finally we have reported in Fig.|Hl the Rayleigh power of the stochastic signal 
from the cosmological background for rjcw = 0.01. For comparison we have 
also reported the noise expected in the VIRGO experiment (Cuoco et al. 1998) 
averaged over one year of integration time. 



6 Discussion 

The two central results of this paper are illustrated by Fig. and Fig. |S1 where 
estimates of the neutrino and gravitational waves diffuse background are re- 
ported, the diffuse neutrino background has been calculated accurately includ- 
ing many details that were not considered in previous calculations. However the 
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final estimate given in this paper is in very reasonable agreement with those of 
previous calculations, indicating that the order of magnitude predictions were 
rather robust. The comparison of this predictions with the best upper limits 
obtained so far by AMANDA for the northern sky, indicates that an increase in 
sensitivity of 2-3 order of magnitude will be required in order to set a stringent 
constraint on the more plausible models of neutrino emission by cosmological 
GRB's. On the other side the detection of extremely powerful bursts of neutri- 
nos in a km'^ detector cannot be ruled out, given the extreme variability of the 
energetic among the various bursts. 

In SjSlit has been shown that the GW emission from single bursts at cosmo- 
logical distances, if the 7-ray prompt emission is beamed with a small angle as 
suggested by afterglow observations, is expected to be c?l /i < 10"^'' Gpc/-\/Hz 
which is well below the detection threshold of presently planned experiments. 
However if the emission is beamed only a small fraction (one over 300) of the 
GRB's are observed by 7-ray satellites. This implies that small amplitude pulse 
trains of GW's impinge over the detector at a frequency that could be as high 
as half per minute, which is practically a continuous signal. Even if the individ- 
ual pulse is small, integrating over a reasonable observation time (order of one 
year) an excess Rayleigh power should emerge from the instrumental noise. The 
frequency spectrum of this eventual excess should give a direct information on 
the characteristic time scale of the collapse and on the cosmological evolution 
of the GRB rate in the recent past {z w 1 — 5). The predicted amplitude is 
conservatively estimated to be in the range of 5 x 10~^^ 1/-\/Hz averaging the 
Rayleigh power over one year. This estimate is rather robust because does not 
depends on the beaming factor and depends only slightly from the large scale 
distribution of the GRB sources. The cosmological origin of the excess noise 
could be proved by detecting a dipole anisotropy. In addition the auto correla- 
tion spectrum of the noise shall carry an imprint of the characteristic duration of 
the GW pulse trains. The detection of those two signatures, even if perhaps not 
possible with presently planned experiments, can give the important additional 
evidence of the cosmological origin of the stochastic signal and informations on 
the physics of GRB's. 
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